(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 45848, 1295]*) (*NotebookOutlinePosition[ 46828, 1325]*) (* CellTagsIndexPosition[ 46784, 1321]*) (*WindowFrame->Normal*) Notebook[{ Cell["\:0412\:044b\:0447\:0438\:0441\:043b\:0438\:043c \:043c\:043e\:043c\ \:0435\:043d\:0442\:044b \:0438\:043d\:0435\:0440\:0446\:0438\:0438", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(TraditionalForm\`Cz = m*r^2/2\), "\[IndentingNewLine]", \(TraditionalForm\`Cxy = m*\((h^2/12 + r^2/4)\)\)}], "Input"], Cell[BoxData[ \(\(m\ r\^2\)\/2\)], "Output"], Cell[BoxData[ \(m\ \((h\^2\/12 + r\^2\/4)\)\)], "Output"] }, Closed]], Cell["\:0413\:0438\:0440\:043e\:0441\:043a\:043e\:043f\:0438\:0447\:0435\:0441\ \:043a\:0438\:0439 \:043c\:043e\:043c\:0435\:043d\:0442", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Hg = Cz*\[Omega]\)], "Input"], Cell[BoxData[ \(1\/2\ m\ r\^2\ \[Omega]\)], "Output"] }, Closed]], Cell["\:041e\:0431\:043e\:0431\:0449\:0435\:043d\:043d\:044b\:0435 \:043a\ \:043e\:043e\:0440\:0434\:0438\:043d\:0430\:0442\:044b", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(TraditionalForm\`s = D[\[Beta][t], t]\), "\[IndentingNewLine]", \(TraditionalForm\`s1 = D[\[Alpha][t], t]\), "\[IndentingNewLine]", \(TraditionalForm\`s2 = D[\[Gamma][t], t]\)}], "Input"], Cell[BoxData[ RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]], "Output"], Cell[BoxData[ RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]], "Output"], Cell[BoxData[ RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]], "Output"] }, Closed]], Cell["\:0421\:0438\:043b\:044b, \:0434\:0435\:0439\:0441\:0442\:0432\:0443\ \:044e\:0449\:0438\:0435 \:043f\:043e \:043e\:0441\:044f\:043c", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(fx = 0\), "\[IndentingNewLine]", \(fy = \((Hg\ s1 - d\ s)\)\ Cos[\[Beta][t]]\), "\[IndentingNewLine]", \(fz = 0\)}], "Input"], Cell[BoxData[ \(0\)], "Output"], Cell[BoxData[ RowBox[{\(Cos[\[Beta][t]]\), " ", RowBox[{"(", RowBox[{ RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"d", " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], ")"}]}]], "Output"], Cell[BoxData[ \(0\)], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(fxyz = {{fx}, {fy}, {fz}} // MatrixForm\)\(\[IndentingNewLine]\) \)\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0"}, { RowBox[{\(Cos[\[Beta][t]]\), " ", RowBox[{"(", RowBox[{ RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"d", " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], ")"}]}]}, {"0"} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Closed]], Cell["\:041a\:0438\:043d\:0435\:0447\:0435\:0441\:043a\:0430\:044f \:044d\ \:043d\:0435\:0440\:0433\:0438\:044f", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(Tk = Cxy*s^2/2 + Cxy*s1^2/2 + Cz*\((\[Omega] + s2)\)^2/2\)], "Input"], Cell[BoxData[ RowBox[{ RowBox[{\(1\/2\), " ", "m", " ", \((h\^2\/12 + r\^2\/4)\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}], "+", RowBox[{\(1\/2\), " ", "m", " ", \((h\^2\/12 + r\^2\/4)\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "2"]}], "+", RowBox[{\(1\/4\), " ", "m", " ", \(r\^2\), " ", SuperscriptBox[ RowBox[{"(", RowBox[{"\[Omega]", "+", RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], ")"}], "2"]}]}]], "Output"] }, Closed]], Cell["\:0417\:0430\:043f\:0438\:0441\:044c \:0443\:0440\:0430\:0432\:043d\ \:0435\:043d\:0438\:0439 \:041b\:0430\:0433\:0440\:0430\:043d\:0436\:0430", \ "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(wx = Integrate\ [fx, \ \[Alpha][t]]\), "\[IndentingNewLine]", \(wy = Integrate\ [fy, \[Beta][t]]\), "\[IndentingNewLine]", \(wz = Integrate[fz, \ \[Gamma][t]]\)}], "Input"], Cell[BoxData[ \(0\)], "Output"], Cell[BoxData[ RowBox[{ RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Sin[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"d", " ", \(Sin[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]], "Output"], Cell[BoxData[ \(0\)], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(qxyz = {{D[wx, \[Alpha][t]]}, {D[wy, \[Beta][t]]}, {D[ wz, \[Gamma][t]]}}\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{\({0}\), ",", RowBox[{"{", RowBox[{ RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{"d", " ", \(Cos[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}], ",", \({0}\)}], "}"}]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(ddTkDqdt = {{D[D[Tk, s1], t]}, {D[D[Tk, s], t]}, {D[D[Tk, s2], t]}}\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}]}], "}"}]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(dTkDq = {{D[Tk, \[Alpha][t]]}, {D[Tk, \[Beta][t]]}, {D[ Tk, \[Gamma][t]]}}\)], "Input"], Cell[BoxData[ \({{0}, {0}, {0}}\)], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(df = ddTkDqdt - dTkDq - qxyz\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{\(-\(1\/2\)\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"d", " ", \(Cos[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}]}], "}"}]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(df = ReplaceAll[df, D[s1, t] \[Rule] D[q11[t], t]]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["q11", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{\(-\(1\/2\)\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"d", " ", \(Cos[\[Beta][t]]\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "+", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}]}], "}"}]], "Output"] }, Closed]], Cell["\:041f\:0440\:0435\:0434\:0441\:0442\:0430\:0432\:043b\:0435\:043d\:0438\ \:0435 \:0432 \:0444\:043e\:0440\:043c\:0435 \:041a\:043e\:0448\:0438", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(m1 = {Solve[q11[t] - D[\[Alpha][t], t] \[Equal] 0, D[\[Alpha][t], t]], \ Solve[Flatten[df[\([1, 1]\)]] \[Equal] 0, D[q11[t], t]]} // MatrixForm\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(q11[t]\)}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]}, { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["q11", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", "0"}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(df = ReplaceAll[ df, \ {D[\[Beta][t], t] \[Rule] q21[t], \ \[Beta][t] \[Rule] q22[t], D[\[Beta][t], {t, 2}] \[Rule] D[q21[t], t]}]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["q11", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{\(d\ Cos[q22[t]]\ q21[t]\), "+", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["q21", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[q22[t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}]}], "}"}]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(m2 = {{q21[t] - D[\[Beta][t], t]}, {D[q22[t], t] - D[\[Beta][t], t]}, \ {Flatten[df[\([2, 1]\)]]}} // MatrixForm\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{\(q21[t]\), "-", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}, { RowBox[{ RowBox[{ SuperscriptBox["q22", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}, { RowBox[{\(d\ Cos[q22[t]]\ q21[t]\), "+", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["q21", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[q22[t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}]} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(m2k = {Solve[q21[t] - D[\[Beta][t], t] \[Equal] 0, D[\[Beta][t], t]], Solve\ [D[q22[t], t] - D[\[Beta][t], t] \[Equal] 0, \ D[q22[t], t]], \ Solve[Flatten[df[\([2, 1]\)]] \[Equal] 0, D[q21[t], t]]} // MatrixForm\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(q21[t]\)}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]}, { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["q22", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]}, { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["q21", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{\(2\ d\ Cos[q22[t]]\ q21[t]\), "-", RowBox[{ "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[q22[t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], ")"}]}], \(m\ \((h\^2 + 3\ r\^2)\)\)]}]}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(df = ReplaceAll[df, \ {D[s2, t] \[Rule] D[q32[t], \ t]}]\)], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["q11", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{\(d\ Cos[q22[t]]\ q21[t]\), "+", RowBox[{"m", " ", \((h\^2\/12 + r\^2\/4)\), " ", RowBox[{ SuperscriptBox["q21", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "-", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", "\[Omega]", " ", \(Cos[q22[t]]\), " ", RowBox[{ SuperscriptBox["\[Alpha]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{\(1\/2\), " ", "m", " ", \(r\^2\), " ", RowBox[{ SuperscriptBox["q32", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}], "}"}]}], "}"}]], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(m3 = {Solve[q31[t] - D[\[Gamma][t], t] \[Equal] 0, D[\[Gamma][t], t]], \ Solve[Flatten[df[\([3, 1]\)]] \[Equal] 0, D[q32[t], t]]} // MatrixForm\)], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["\[Gamma]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", \(q31[t]\)}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]}, { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["q32", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", "0"}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Closed]], Cell["\:0412\:0432\:043e\:0434 \:0443\:043f\:0440\:043e\:0449\:0435\:043d\ \:0438\:0439", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(ReplaceAll[m1, D[\[Alpha][t], t] \[Rule] 0]\), "\[IndentingNewLine]", \(ReplaceAll[m3, D[\[Gamma][t], t] \[Rule] 0]\), "\[IndentingNewLine]", \(m2 = ReplaceAll[ m2, \ {Cos[q22[t]] \[Rule] 1, \ D[q21[t], t] \[Rule] 0, D[\[Alpha][t], t] \[Rule] \[CapitalOmega]}]\)}], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ {\(0 \[Rule] q11[t]\)} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]}, { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["q11", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", "0"}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ {\(0 \[Rule] q31[t]\)} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]}, { RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{ RowBox[{ SuperscriptBox["q32", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "\[Rule]", "0"}]} }, RowSpacings->1, ColumnAlignments->{Left}], Column], "\[NoBreak]", ")"}]} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{\(q21[t]\), "-", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}, { RowBox[{ RowBox[{ SuperscriptBox["q22", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}, {\(\(-\(1\/2\)\)\ m\ r\^2\ \[Omega]\ \[CapitalOmega] + d\ q21[t]\)} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"] }, Closed]], Cell["\:0420\:0435\:0448\:0435\:043d\:0438\:0435", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(resh0 = Flatten[Solve[Flatten[m2[\([1, 3]\)]] \[Equal] 0, q21[t]]]\), "\[IndentingNewLine]", \(Clear[\[Beta]]\), "\[IndentingNewLine]", \(m2 = ReplaceAll[m2, resh0[\([1]\)]]\), "\[IndentingNewLine]", \(resh = Flatten[DSolve[{\(m2[\([1, 1]\)]\)[\([1]\)] \[Equal] 0, \[Beta][0] \[Equal] 0}, \[Beta][t], t]]\)}], "Input"], Cell[BoxData[ \({q21[ t] \[Rule] \(m\ r\^2\ \[Omega]\ \[CapitalOmega]\)\/\(2\ d\)}\)], \ "Output"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{\(\(m\ r\^2\ \[Omega]\ \[CapitalOmega]\)\/\(2\ d\)\), "-", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}, { RowBox[{ RowBox[{ SuperscriptBox["q22", "\[Prime]", MultilineFunction->None], "[", "t", "]"}], "-", RowBox[{ SuperscriptBox["\[Beta]", "\[Prime]", MultilineFunction->None], "[", "t", "]"}]}]}, {"0"} }, RowSpacings->1, ColumnSpacings->1, ColumnAlignments->{Left}], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Output"], Cell[BoxData[ \({\[Beta][ t] \[Rule] \(m\ r\^2\ t\ \[Omega]\ \[CapitalOmega]\)\/\(2\ d\)}\)], \ "Output"] }, Closed]], Cell["\:041f\:043e\:0434\:0441\:0442\:0430\:043d\:043e\:0432\:043a\:0430", \ "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(resh = ReplaceAll[ resh, {m -> dr^2 Pi/4*h*ro, r^2 \[Rule] dr^2/4}]\), "\[IndentingNewLine]", \(resh = ReplaceAll[ resh, {dr \[Rule] 2, \ h \[Rule] 1, d \[Rule] 10000, \[Omega] \[Rule] 20000\ 2\ Pi/60, ro \[Rule] 7.8}]\), "\[IndentingNewLine]", \(resh = ReplaceAll[resh, \[CapitalOmega] \[Rule] 0.2]\)}], "Input"], Cell[BoxData[ \({\[Beta][ t] \[Rule] \(dr\^4\ h\ \[Pi]\ ro\ t\ \[Omega]\ \ \[CapitalOmega]\)\/\(32\ d\)}\)], "Output"], Cell[BoxData[ \({\[Beta][t] \[Rule] 2.5660971442832325`\ t\ \[CapitalOmega]}\)], "Output"], Cell[BoxData[ \({\[Beta][t] \[Rule] 0.5132194288566465`\ t}\)], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(\:041f\:043e\:0441\:0442\:0440\:043e\:0435\:043d\:0438\:0435\ \:0433\ \:0440\:0430\:0444\:0444\:0438\:043a\:0430\)], "Input"], Cell[BoxData[ \(\:0433\:0440\:0430\:0444\:0444\:0438\:043a\:0430\ \:041f\:043e\:0441\ \:0442\:0440\:043e\:0435\:043d\:0438\:0435\)], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[{ \(tn = 0\), "\[IndentingNewLine]", \(tk = 5\)}], "Input"], Cell[BoxData[ \(0\)], "Output"], Cell[BoxData[ \(5\)], "Output"] }, Closed]], Cell[CellGroupData[{ Cell[BoxData[ \(Plot[resh[\([1, 2]\)], \ {t, \ tn, \ tk}]\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.190476 0.0147151 0.229377 [ [.21429 .00222 -3 -9 ] [.21429 .00222 3 0 ] [.40476 .00222 -3 -9 ] [.40476 .00222 3 0 ] [.59524 .00222 -3 -9 ] [.59524 .00222 3 0 ] [.78571 .00222 -3 -9 ] [.78571 .00222 3 0 ] [.97619 .00222 -3 -9 ] [.97619 .00222 3 0 ] [.01131 .1294 -18 -4.5 ] [.01131 .1294 0 4.5 ] [.01131 .24409 -6 -4.5 ] [.01131 .24409 0 4.5 ] [.01131 .35878 -18 -4.5 ] [.01131 .35878 0 4.5 ] [.01131 .47347 -6 -4.5 ] [.01131 .47347 0 4.5 ] [.01131 .58816 -18 -4.5 ] [.01131 .58816 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .21429 .01472 m .21429 .02097 L s [(1)] .21429 .00222 0 1 Mshowa .40476 .01472 m .40476 .02097 L s [(2)] .40476 .00222 0 1 Mshowa .59524 .01472 m .59524 .02097 L s [(3)] .59524 .00222 0 1 Mshowa .78571 .01472 m .78571 .02097 L s [(4)] .78571 .00222 0 1 Mshowa .97619 .01472 m .97619 .02097 L s [(5)] .97619 .00222 0 1 Mshowa .125 Mabswid .0619 .01472 m .0619 .01847 L s .1 .01472 m .1 .01847 L s .1381 .01472 m .1381 .01847 L s .17619 .01472 m .17619 .01847 L s .25238 .01472 m .25238 .01847 L s .29048 .01472 m .29048 .01847 L s .32857 .01472 m .32857 .01847 L s .36667 .01472 m .36667 .01847 L s .44286 .01472 m .44286 .01847 L s .48095 .01472 m .48095 .01847 L s .51905 .01472 m .51905 .01847 L s .55714 .01472 m .55714 .01847 L s .63333 .01472 m .63333 .01847 L s .67143 .01472 m .67143 .01847 L s .70952 .01472 m .70952 .01847 L s .74762 .01472 m .74762 .01847 L s .82381 .01472 m .82381 .01847 L s .8619 .01472 m .8619 .01847 L s .9 .01472 m .9 .01847 L s .9381 .01472 m .9381 .01847 L s .25 Mabswid 0 .01472 m 1 .01472 L s .02381 .1294 m .03006 .1294 L s [(0.5)] .01131 .1294 1 0 Mshowa .02381 .24409 m .03006 .24409 L s [(1)] .01131 .24409 1 0 Mshowa .02381 .35878 m .03006 .35878 L s [(1.5)] .01131 .35878 1 0 Mshowa .02381 .47347 m .03006 .47347 L s [(2)] .01131 .47347 1 0 Mshowa .02381 .58816 m .03006 .58816 L s [(2.5)] .01131 .58816 1 0 Mshowa .125 Mabswid .02381 .03765 m .02756 .03765 L s .02381 .06059 m .02756 .06059 L s .02381 .08353 m .02756 .08353 L s .02381 .10647 m .02756 .10647 L s .02381 .15234 m .02756 .15234 L s .02381 .17528 m .02756 .17528 L s .02381 .19822 m .02756 .19822 L s .02381 .22115 m .02756 .22115 L s .02381 .26703 m .02756 .26703 L s .02381 .28997 m .02756 .28997 L s .02381 .31291 m .02756 .31291 L s .02381 .33584 m .02756 .33584 L s .02381 .38172 m .02756 .38172 L s .02381 .40466 m .02756 .40466 L s .02381 .42759 m .02756 .42759 L s .02381 .45053 m .02756 .45053 L s .02381 .49641 m .02756 .49641 L s .02381 .51934 m .02756 .51934 L s .02381 .54228 m .02756 .54228 L s .02381 .56522 m .02756 .56522 L s .02381 .6111 m .02756 .6111 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .01472 m .06244 .03859 L .10458 .06463 L .14415 .08909 L .18221 .11261 L .22272 .13765 L .26171 .16175 L .30316 .18736 L .34309 .21204 L .3815 .23578 L .42237 .26104 L .46172 .28536 L .49955 .30874 L .53984 .33364 L .57861 .3576 L .61984 .38308 L .65954 .40762 L .69774 .43123 L .73838 .45635 L .77751 .48053 L .81909 .50623 L .85916 .53099 L .89771 .55482 L .93871 .58016 L .97619 .60332 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgP0000000000001V]S40oooo00@0fi0j0000000003X0T=_o;`3oool0 103KT3X000000000000jT=/60?ooo`00B03oool00`3o]VH0>Y3K0?ooo`0`0?ooo`03090j000006H0 ]_oo02l0oooo00H0om^@03XjT03oool0oonf06H0>P2@fola0?ooo`030?nfIP0jT=/0oooo02h0oooo 00H0om^@03XjT03oool0oooo090j>P2@fol50?ooo`00B03oool00`3o]VH0>Y3K0?ooo`0a0?ooo`03 0?nfIP000000I[Ko0300oooo00<0]VH00000001V]_l0;`3oool01P3ookH0IP00000000000000000j 093Koc80oooo00<0]VHj093Ko`3oool0103oool004P0oooo00<0okIV03Z@f`3oool0Y00oooo02l0oooo00@0]VH006Jfo`3o]VH0>Y3K<03oool0 102fIP00000000000000I[H60?ooo`00A`3oool00`2@>VH0][IV03Z@f`0`0?ooo`060?oKT00j>Y00 oooo0?oo]P1V06H0]_oo;@3oool01P3ookH0IP0j093Ko`3ookH0IP1V0;Koobl0oooo00@0oonf06H0 >P2@T6H0>Y3K<03oool00`2fISX0T=_o0?ooo`070?ooo`00A`3oool00`3o]VH0000003Z@f`0a0?oo o`040=^@>P000000000003Z@fbl0oooo00@0om^@03X000000000I[Ko<@3oool00`3ofi00>P0003Z@ f`0`0?ooo`040;IV000000000000001V]PH0oooo003o0?ooob40oooo000E0?ooo`030000003oool0 oooo0?l0oooo2@3oool001D0oooo00<000000?ooo`3oool0o`3oool90?ooo`005@3oool00`000000 oooo0?ooo`3o0?ooo`T0oooo000E0?ooo`030000003oool0oooo0?l0oooo2@3oool000h0ooooo`00 000C000000005@3oool3000000L0oooo00<000000?ooo`3oool0203oool00`000000oooo0?ooo`07 0?ooo`030000003oool0oooo00L0oooo00<000000?ooo`3oool0203oool00`000000oooo0?ooo`07 0?ooo`030000003oool0oooo00L0oooo00<000000?ooo`3oool0203oool00`000000oooo0?ooo`07 0?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool01`3oool00`000000oooo0?ooo`07 0?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool01`3oool00`000000oooo0?ooo`07 0?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool01`3oool00`000000oooo0?ooo`07 0?ooo`030000003oool0oooo00P0oooo00<000000?ooo`3oool01`3oool00`000000oooo0?ooo`08 0?ooo`030000003oool0oooo00L0oooo00<000000?ooo`3oool01`3oool00`000000oooo0?ooo`08 0?ooo`030000003oool0oooo00L0oooo00<000000?ooo`3oool01@3oool001D0oooo00@000000?oo o`3oool00000<03oool00`000000oooo0?ooo`0a0?ooo`030000003oool0oooo0300oooo00<00000 0?ooo`3oool0<@3oool00`000000oooo0?ooo`0a0?ooo`030000003oool0oooo00D0oooo000E0?oo o`050000003oool0oooo0?ooo`000000o`3oool70?ooo`005@3oool00`000000oooo0?ooo`020?oo o`800000o`3oool50?ooo`005@3oool00`000000oooo0?ooo`040?ooo`030000003oool0oooo0?l0 oooo0P3oool001D0oooo0P0000060?ooo`800000o`3oool20?ooo`005@3oool00`000000oooo0?oo o`070?ooo`030000003oool0oooo0?h0oooo000E0?ooo`030000003oool0oooo00P0oooo0P00003n 0?ooo`005@3oool00`000000oooo0?ooo`0:0?ooo`800000o03oool001D0oooo00<000000?ooo`3o ool0303oool00`000000oooo0?ooo`3i0?ooo`005@3oool00`000000oooo0?ooo`0=0?ooo`800000 n@3oool001D0oooo00<000000?ooo`3oool03`3oool200000?L0oooo000E0?ooo`8000004P3oool2 00000?D0oooo000E0?ooo`030000003oool0oooo01<0oooo00<000000?ooo`3oool0lP3oool001D0 oooo00<000000?ooo`3oool0503oool200000?80oooo000E0?ooo`030000003oool0oooo01H0oooo 0P00003`0?ooo`005@3oool00`000000oooo0?ooo`0H0?ooo`800000kP3oool001D0oooo00<00000 0?ooo`3oool06P3oool200000>`0oooo000E0?ooo`8000007@3oool200000>X0oooo000E0?ooo`03 0000003oool0oooo01h0oooo00<000000?ooo`3oool0i`3oool001D0oooo00<000000?ooo`3oool0 7`3oool200000>L0oooo000E0?ooo`030000003oool0oooo0240oooo00<000000?ooo`3oool0i03o ool001D0oooo00<000000?ooo`3oool08P3oool00`000000oooo0?ooo`3S0?ooo`005@3oool00`00 0000oooo0?ooo`0S0?ooo`800000h`3oool001D0oooo0P00000V0?ooo`030000003oool0oooo0>00 oooo000E0?ooo`030000003oool0oooo02H0oooo0P00003P0?ooo`005@3oool00`000000oooo0?oo o`0X0?ooo`030000003oool0oooo0=d0oooo000E0?ooo`030000003oool0oooo02T0oooo0P00003M 0?ooo`0000D0oooo0=^@>P000000000003Z@f`020?ooo`030?oo]P1V000003Z@00<0oooo00@0fi0j 000000000000>Y3K103oool00`000000oooo0?ooo`0[0?ooo`030000003oool0oooo0=X0oooo0000 1P3oool0T6J@0?ooo`3oool0T3YV0;Koo`H0oooo00H0om^@03XjT03oool0oooo090j>P2@fol30?oo o`030000003oool0oooo02`0oooo0P00003J0?ooo`0000H0oonf06HjT03oool0oooo0;IV>P2@fol: 0?ooo`030;IV>P2@fol0oooo0080oooo0`00000^0?ooo`800000f03oool000060?oo]P1V>Y00oooo 0?ooo`2fISX0T=_o1`3oool0102fIP00000000000000I[H40?ooo`030000003oool0oooo0300oooo 00<000000?ooo`3oool0e@3oool000060?ooo`2@>SX0T=_o0?oo]P1V06H0]_oo1`3oool00`2fISX0 T=_o0?ooo`050?ooo`030000003oool0oooo0340oooo0P00003E0?ooo`0000D0oooo0?oKT00j0000 000006Jfo`080?ooo`040;IV000000000000001V]P@0oooo00<000000?ooo`3oool0<`3oool00`00 0000oooo0?ooo`3B0?ooo`005@3oool00`000000oooo0?ooo`0d0?ooo`800000dP3oool001D0oooo 00<000000?ooo`3oool0=P3oool200000=00oooo000E0?ooo`030000003oool0oooo03P0oooo0P00 003>0?ooo`005@3oool2000003/0oooo0P00003<0?ooo`005@3oool00`000000oooo0?ooo`0l0?oo o`800000bP3oool001D0oooo00<000000?ooo`3oool0?P3oool00`000000oooo0?ooo`370?ooo`00 5@3oool00`000000oooo0?ooo`0o0?ooo`800000a`3oool001D0oooo00<000000?ooo`3oool0@@3o ool00`000000oooo0?ooo`340?ooo`005@3oool00`000000oooo0?ooo`120?ooo`800000a03oool0 01D0oooo0P0000150?ooo`800000`P3oool001D0oooo00<000000?ooo`3oool0AP3oool00`000000 oooo0?ooo`2o0?ooo`005@3oool00`000000oooo0?ooo`170?ooo`800000_`3oool001D0oooo00<0 00000?ooo`3oool0B@3oool00`000000oooo0?ooo`2l0?ooo`005@3oool00`000000oooo0?ooo`1: 0?ooo`800000_03oool001D0oooo00<000000?ooo`3oool0C03oool00`000000oooo0?ooo`2i0?oo o`005@3oool2000004h0oooo0P00002i0?ooo`005@3oool00`000000oooo0?ooo`1?0?ooo`800000 ]`3oool001D0oooo00<000000?ooo`3oool0D@3oool00`000000oooo0?ooo`2d0?ooo`005@3oool0 0`000000oooo0?ooo`1B0?ooo`800000]03oool001D0oooo00<000000?ooo`3oool0E03oool00`00 0000oooo0?ooo`2a0?ooo`005@3oool00`000000oooo0?ooo`1E0?ooo`800000/@3oool001D0oooo 0P00001H0?ooo`800000[`3oool001D0oooo00<000000?ooo`3oool0F@3oool00`000000oooo0?oo o`2/0?ooo`005@3oool00`000000oooo0?ooo`1J0?ooo`800000[03oool001D0oooo00<000000?oo o`3oool0G03oool200000:X0oooo000E0?ooo`030000003oool0oooo05h0oooo00<000000?ooo`3o ool0Y`3oool000`0oooo00H0oonf06H00000000000000000>P2@fol30?ooo`030000003oool0oooo 05l0oooo0P00002W0?ooo`003P3oool00`3o]VH0>Y3K0?ooo`040?ooo`030000003oool0oooo0640 oooo00<000000?ooo`3oool0Y03oool000h0oooo00<0okIV03Z@f`3oool0103oool300000680oooo 0P00002T0?ooo`003P3oool00`3o]VH0>Y3K0?ooo`040?ooo`030000003oool0oooo06@0oooo0P00 002R0?ooo`003@3oool00`2@>VH0][IV03Z@f`050?ooo`030000003oool0oooo06H0oooo00<00000 0?ooo`3oool0W`3oool000d0oooo00<0okIV0000000jT=/01@3oool00`000000oooo0?ooo`1W0?oo o`800000W`3oool001D0oooo00<000000?ooo`3oool0J@3oool00`000000oooo0?ooo`2L0?ooo`00 5@3oool00`000000oooo0?ooo`1Z0?ooo`800000W03oool001D0oooo0P00001]0?ooo`030000003o ool0oooo09T0oooo000E0?ooo`030000003oool0oooo06d0oooo0P00002I0?ooo`005@3oool00`00 0000oooo0?ooo`1_0?ooo`800000U`3oool001D0oooo00<000000?ooo`3oool0L@3oool00`000000 oooo0?ooo`2D0?ooo`005@3oool00`000000oooo0?ooo`1b0?ooo`800000U03oool001D0oooo00<0 00000?ooo`3oool0M03oool00`000000oooo0?ooo`2A0?ooo`005@3oool2000007H0oooo0P00002A 0?ooo`005@3oool00`000000oooo0?ooo`1g0?ooo`800000S`3oool001D0oooo00<000000?ooo`3o ool0N@3oool00`000000oooo0?ooo`2<0?ooo`005@3oool00`000000oooo0?ooo`1j0?ooo`800000 S03oool001D0oooo00<000000?ooo`3oool0O03oool2000008X0oooo000E0?ooo`030000003oool0 oooo07h0oooo00<000000?ooo`3oool0Q`3oool001D0oooo0P0000200?ooo`800000Q`3oool001D0 oooo00<000000?ooo`3oool0P@3oool00`000000oooo0?ooo`240?ooo`005@3oool00`000000oooo 0?ooo`220?ooo`800000Q03oool001D0oooo00<000000?ooo`3oool0Q03oool200000880oooo000E 0?ooo`030000003oool0oooo08H0oooo00<000000?ooo`3oool0O`3oool001D0oooo00<000000?oo o`3oool0Q`3oool2000007l0oooo000E0?ooo`030000003oool0oooo08T0oooo00<000000?ooo`3o ool0O03oool001D0oooo0P00002;0?ooo`800000O03oool001D0oooo00<000000?ooo`3oool0S03o ool2000007X0oooo000E0?ooo`030000003oool0oooo08h0oooo0P00001h0?ooo`005@3oool00`00 0000oooo0?ooo`2@0?ooo`800000MP3oool0000:0?oo]P1V000000000000000003X0T=_o0?ooo`3o okH0IP00000jT0<0oooo00@0fi0j000000000000>Y3K103oool00`000000oooo0?ooo`2B0?ooo`80 0000M03oool00080oooo00<0okIV03Z@f`3oool01`3oool01P3ofi00>SZ@0?ooo`3oool0T3Xj093K o`<0oooo00<000000?ooo`3oool0U03oool00`000000oooo0?ooo`1a0?ooo`000P3oool00`3o]VH0 >Y3K0?ooo`0;0?ooo`030;IV>P2@fol0oooo0080oooo0`00002E0?ooo`800000L@3oool00080oooo 00<0okIV03Z@f`3oool0203oool0102fIP00000000000000I[H40?ooo`030000003oool0oooo09L0 oooo00<000000?ooo`3oool0KP3oool000040?ooo`2@>VH0][IV03Z@f`T0oooo00<0]VHj093Ko`3o ool01@3oool00`000000oooo0?ooo`2H0?ooo`800000KP3oool000040?ooo`3o]VH0000003Z@f`T0 oooo00@0]VH000000000000006Jf103oool00`000000oooo0?ooo`2J0?ooo`800000K03oool001D0 oooo00<000000?ooo`3oool0W03oool00`000000oooo0?ooo`1Y0?ooo`005@3oool00`000000oooo 0?ooo`2M0?ooo`800000J@3oool001D0oooo0P00002P0?ooo`030000003oool0oooo06H0oooo000E 0?ooo`030000003oool0oooo0:00oooo0P00001V0?ooo`005@3oool00`000000oooo0?ooo`2R0?oo o`030000003oool0oooo06<0oooo000E0?ooo`030000003oool0oooo0:<0oooo0P00001S0?ooo`00 5@3oool00`000000oooo0?ooo`2U0?ooo`800000H@3oool001D0oooo00<000000?ooo`3oool0Y`3o ool00`000000oooo0?ooo`1N0?ooo`005@3oool200000:T0oooo0P00001N0?ooo`005@3oool00`00 0000oooo0?ooo`2Z0?ooo`030000003oool0oooo05/0oooo000E0?ooo`030000003oool0oooo0:/0 oooo0P00001K0?ooo`005@3oool00`000000oooo0?ooo`2]0?ooo`800000F@3oool001D0oooo00<0 00000?ooo`3oool0[`3oool00`000000oooo0?ooo`1F0?ooo`005@3oool00`000000oooo0?ooo`2` 0?ooo`800000EP3oool001D0oooo0P00002c0?ooo`800000E03oool001D0oooo00<000000?ooo`3o ool0]03oool00`000000oooo0?ooo`1A0?ooo`005@3oool00`000000oooo0?ooo`2e0?ooo`800000 D@3oool001D0oooo00<000000?ooo`3oool0]`3oool00`000000oooo0?ooo`1>0?ooo`005@3oool0 0`000000oooo0?ooo`2h0?ooo`800000CP3oool001D0oooo00<000000?ooo`3oool0^P3oool20000 04`0oooo000E0?ooo`030000003oool0oooo0;`0oooo00<000000?ooo`3oool0B@3oool001D0oooo 0P00002n0?ooo`800000B@3oool001D0oooo00<000000?ooo`3oool0_`3oool00`000000oooo0?oo o`160?ooo`005@3oool00`000000oooo0?ooo`300?ooo`800000AP3oool001D0oooo00<000000?oo o`3oool0`P3oool00`000000oooo0?ooo`130?ooo`00303oool01P3o]VH00000000000000000001V 0;Koo`<0oooo00<000000?ooo`3oool0``3oool2000004<0oooo000=0?ooo`03090j000006H0]_oo 00D0oooo00<000000?ooo`3oool0a@3oool200000440oooo000>0?ooo`030?nfIP000000I[Ko00@0 oooo0`0000370?ooo`030000003oool0oooo03h0oooo000?0?ooo`030?oo]P2f]_l0oooo00<0oooo 00<000000?ooo`3oool0b03oool2000003h0oooo000<0?ooo`060?oKT00j>Y00oooo0?oo]P1V06H0 ]_oo0`3oool00`000000oooo0?ooo`3:0?ooo`030000003oool0oooo03/0oooo000=0?ooo`040=^@ >P000000000003Z@f`@0oooo00<000000?ooo`3oool0b`3oool2000003/0oooo000E0?ooo`030000 003oool0oooo003oool001D0oooo00<000000?ooo`3oool0cP3o ool2000003P0oooo000E0?ooo`800000d@3oool2000003H0oooo000E0?ooo`030000003oool0oooo 0=80oooo00<000000?ooo`3oool0<`3oool001D0oooo00<000000?ooo`3oool0d`3oool2000003<0 oooo000E0?ooo`030000003oool0oooo0=D0oooo00<000000?ooo`3oool0<03oool001D0oooo00<0 00000?ooo`3oool0eP3oool200000300oooo000E0?ooo`030000003oool0oooo0=P0oooo00<00000 0?ooo`3oool0;@3oool001D0oooo0P00003J0?ooo`800000;@3oool001D0oooo00<000000?ooo`3o ool0f`3oool2000002/0oooo000E0?ooo`030000003oool0oooo0=d0oooo00<000000?ooo`3oool0 :03oool001D0oooo00<000000?ooo`3oool0gP3oool2000002P0oooo000E0?ooo`030000003oool0 oooo0>00oooo00<000000?ooo`3oool09@3oool001D0oooo00<000000?ooo`3oool0h@3oool20000 02D0oooo000E0?ooo`030000003oool0oooo0><0oooo0P00000S0?ooo`005@3oool200000>H0oooo 0P00000Q0?ooo`005@3oool00`000000oooo0?ooo`3W0?ooo`8000007`3oool001D0oooo00<00000 0?ooo`3oool0j@3oool2000001d0oooo000E0?ooo`030000003oool0oooo0>/0oooo00<000000?oo o`3oool06P3oool001D0oooo00<000000?ooo`3oool0k03oool2000001X0oooo000E0?ooo`030000 003oool0oooo0>h0oooo00<000000?ooo`3oool05`3oool001D0oooo0P00003`0?ooo`8000005`3o ool001D0oooo00<000000?ooo`3oool0l@3oool2000001D0oooo000E0?ooo`030000003oool0oooo 0?<0oooo00<000000?ooo`3oool04P3oool001D0oooo00<000000?ooo`3oool0m03oool200000180 oooo00002P3o]VH00000000000000000001V0;Koo`3oool0oonf06H00000>Y030?ooo`040=^@>P00 0000000003Z@f`@0oooo00<000000?ooo`3oool0mP3oool00`000000oooo0?ooo`0?0?ooo`0000@0 oooo090j000006H0]_oo203oool01P3ofi00>SZ@0?ooo`3oool0T3Xj093Ko`<0oooo00<000000?oo o`3oool0m`3oool2000000l0oooo00020?ooo`030?nfIP000000I[Ko00/0oooo00<0]VHj093Ko`3o ool00P3oool300000?T0oooo0P00000=0?ooo`000`3oool00`3ookH0][Ko0?ooo`070?ooo`040;IV 000000000000001V]P@0oooo00<000000?ooo`3oool0n`3oool00`000000oooo0?ooo`0:0?ooo`00 00H0om^@03XjT03oool0oonf06H0IP2fool70?ooo`030;IV>P2@fol0oooo00D0oooo00<000000?oo o`3oool0o03oool2000000X0oooo00001@3oool0fi0j000000000000>Y3K00P0oooo00@0]VH00000 0000000006Jf103oool00`000000oooo0?ooo`3n0?ooo`800000203oool001D0oooo00<000000?oo o`3oool0o`3oool10?ooo`030000003oool0oooo00D0oooo000E0?ooo`030000003oool0oooo0?l0 oooo2@3oool001D0oooo0P00003o0?ooo`X0oooo000E0?ooo`030000003oool0oooo0?l0oooo2@3o ool001D0oooo00<000000?ooo`3oool0o`3oool90?ooo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.41298, -0.201822, \ 0.0192962, 0.0160237}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(ReplaceAll[resh, {t \[Rule] tk}]\)], "Input"], Cell[BoxData[ \({\[Beta][5] \[Rule] 2.5660971442832325`}\)], "Output"] }, Closed]] }, FrontEndVersion->"5.2 for Microsoft Windows", ScreenRectangle->{{0, 1280}, {0, 927}}, WindowSize->{924, 740}, WindowMargins->{{162, Automatic}, {Automatic, 22}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PageHeaders->{{Inherited, Inherited, Inherited}, {None, Inherited, None}}, PrintingOptions->{"PrintingMargins"->{{28.25, 28.25}, {28.25, 28.25}}, "PrintCellBrackets"->False, "PrintRegistrationMarks"->True, "PrintMultipleHorizontalPages"->False, "FirstPageHeader"->False} ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1754, 51, 152, 1, 33, "Text"], Cell[CellGroupData[{ Cell[1931, 56, 143, 2, 48, "Input"], Cell[2077, 60, 48, 1, 45, "Output"], Cell[2128, 63, 61, 1, 45, "Output"] }, Closed]], Cell[2204, 67, 145, 1, 30, "Text"], Cell[CellGroupData[{ Cell[2374, 72, 49, 1, 30, "Input"], Cell[2426, 75, 57, 1, 42, "Output"] }, Closed]], Cell[2498, 79, 139, 1, 30, "Text"], Cell[CellGroupData[{ Cell[2662, 84, 218, 3, 67, "Input"], Cell[2883, 89, 132, 3, 29, "Output"], Cell[3018, 94, 133, 3, 29, "Output"], Cell[3154, 99, 133, 3, 29, "Output"] }, Closed]], Cell[3302, 105, 148, 1, 30, "Text"], Cell[CellGroupData[{ Cell[3475, 110, 156, 3, 70, "Input"], Cell[3634, 115, 35, 1, 29, "Output"], Cell[3672, 118, 500, 12, 42, "Output"], Cell[4175, 132, 35, 1, 29, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[4247, 138, 108, 2, 47, "Input"], Cell[4358, 142, 963, 25, 79, "Output"] }, Closed]], Cell[5336, 170, 121, 1, 30, "Text"], Cell[CellGroupData[{ Cell[5482, 175, 89, 1, 30, "Input"], Cell[5574, 178, 818, 19, 45, "Output"] }, Closed]], Cell[6407, 200, 160, 2, 30, "Text"], Cell[CellGroupData[{ Cell[6592, 206, 203, 3, 70, "Input"], Cell[6798, 211, 35, 1, 29, "Output"], Cell[6836, 214, 448, 10, 42, "Output"], Cell[7287, 226, 35, 1, 29, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[7359, 232, 118, 2, 27, "Input"], Cell[7480, 236, 629, 15, 42, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[8146, 256, 113, 2, 27, "Input"], Cell[8262, 260, 779, 18, 45, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[9078, 283, 119, 2, 27, "Input"], Cell[9200, 287, 49, 1, 29, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[9286, 293, 61, 1, 27, "Input"], Cell[9350, 296, 1295, 29, 45, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[10682, 330, 83, 1, 27, "Input"], Cell[10768, 333, 1282, 29, 45, "Output"] }, Closed]], Cell[12065, 365, 159, 1, 30, "Text"], Cell[CellGroupData[{ Cell[12249, 370, 209, 4, 30, "Input"], Cell[12461, 376, 1318, 34, 57, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[13816, 415, 189, 4, 27, "Input"], Cell[14008, 421, 1091, 25, 45, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[15136, 451, 169, 3, 27, "Input"], Cell[15308, 456, 1387, 33, 83, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[16732, 494, 288, 5, 47, "Input"], Cell[17023, 501, 2925, 67, 89, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[19985, 573, 89, 1, 27, "Input"], Cell[20077, 576, 1078, 25, 45, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[21192, 606, 209, 4, 27, "Input"], Cell[21404, 612, 1318, 34, 57, "Output"] }, Closed]], Cell[22737, 649, 97, 1, 30, "Text"], Cell[CellGroupData[{ Cell[22859, 654, 331, 6, 70, "Input"], Cell[23193, 662, 1079, 29, 57, "Output"], Cell[24275, 693, 1079, 29, 57, "Output"], Cell[25357, 724, 899, 23, 79, "Output"] }, Closed]], Cell[26271, 750, 58, 0, 30, "Text"], Cell[CellGroupData[{ Cell[26354, 754, 394, 8, 90, "Input"], Cell[26751, 764, 112, 3, 46, "Output"], Cell[26866, 769, 894, 24, 83, "Output"], Cell[27763, 795, 119, 3, 46, "Output"] }, Closed]], Cell[27897, 801, 84, 1, 30, "Text"], Cell[CellGroupData[{ Cell[28006, 806, 408, 10, 70, "Input"], Cell[28417, 818, 132, 3, 46, "Output"], Cell[28552, 823, 103, 2, 29, "Output"], Cell[28658, 827, 77, 1, 29, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[28772, 833, 145, 2, 27, "Input"], Cell[28920, 837, 146, 2, 29, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[29103, 844, 81, 2, 47, "Input"], Cell[29187, 848, 35, 1, 29, "Output"], Cell[29225, 851, 35, 1, 29, "Output"] }, Closed]], Cell[CellGroupData[{ Cell[29297, 857, 74, 1, 27, "Input"], Cell[29374, 860, 16146, 418, 186, 3624, 259, "GraphicsData", "PostScript", \ "Graphics"], Cell[45523, 1280, 130, 3, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[45690, 1288, 65, 1, 30, "Input"], Cell[45758, 1291, 74, 1, 29, "Output"] }, Closed]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)